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PREFACE 


A 1-day set of briefings on the feasibility of assessing crop 
condition and yield from Landsat data was given at NASA Head- 
quarters on September 27, 1977, to allow an assessment of the 
technical status and remaining technical issues on this important 
topic. It allowed an update on the Landsat aspect of yield from 
the earlier 1974 NASA JSC Wheat-Yield Conference (NASA TH 
X-58158, JSC-09256, April 1975). 
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SECTION 1 

INTRODUCTION AND SUMMARY 
By 

Jon D. Erickson and Richard G. Stuff 


The purpose of the briefings was to: 

e Present the technical status of key Investigations In the 
assessment of crop condition and yield using Landsat data. 

e Identify the technical Issues that are currently limiting 
the research progress or applications of Landsat to yield 
estimation. 

e Provide briefing charts and sufficient textural material 
to publish a readable report documenting the technical 
status and Issues. 


The briefing agenda consisted of the same order of subjects 
and speakers as given in the table of contents. 
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The objective of yield Modeling for crop production estimation Is to derive 

a Means of predicting the wlthin-a-year yield and the year-to-year variability 

of yield over some fixed or randonly located unit of area. Yield prediction 

# 

models have traditionally been empirical functions of weather variables (1) or 
In-season sampling of crop dry matter and stand parameters (2). The need for 
Improved yield models Incorporating satellite data was described for the Large 
Area Crop Inventory Experiment (2-4)* and by USOA personnel attending the 
briefings. In addition to better yield predictions for their component role 
In production, they can also contribute to crop Identification and area 
determination by remote sensing since expected values of yield Indicate the 
condition of a standing crop and the probability that an area or portion 
thereof will not be harvested can be computed. Preliminary studies Indicate 
that the r r^ui remen ts for interpreting Landsat data for yield may be suffi- 
ciently similar to those of signature extension that It Is feasible to 
investigate the automated estimation of production (4-24). 

The model approaches proposed for estimating yield from Landsat data are 
based on the explicit or implicit use of crop condition variables. Although 
crop condition and Landsat data may be analyzed separately, it should be 
recognized that the ultimate function is to employ the results in yield or 
production prediction. A preferable method for describing or quantifying 
crop condition is in terms of expected yield per unit area and conversely, 
expected yield relative to normal should provide the best available quanti- 
fication of crop condition. 

To date, none of the crop condition indicators described in the briefings 
have been functionally related to yield in a tested model. Some of the crop 
condition indicators used in the described Landsat-yield studies are detrac- 
tant cause (3-16), percent green cover (4-11), redness in Landsat color 
composites (5-4), threshold index of transformed Landsat data (5-7), stand 
quality (5-28), two or three classes of stress estimated at 12.5 x 12.5 n.m. 
coordinates on Landsat images (7-2), and LAI at given growth stages (6-7). 

A correlation between percent green cover and yields was indicated for 


*Number-numbers in parentheses refer to pages in this report. 
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observations fro* a LACIE test site (4-11). Also the correlations between 
leaf area duration and yield as reported by Wei bank et al. are frequently 
cited (3*1) *ven though these correlations were Insignificantly low or 
negative for winter wheat. In other cases the field measured crop 
.condition-yield comparison Is omitted (6-7). The lack of yield models 
; based on ground observed crop conditiori variables could be considered one 
of the voids in the technology for estimating yields from Landsat. 


The main potential or feasibility indicators for using Landsat data to esti- 
mate yield are the spot correlations (not models) between KSS data and crop 
c-’ditlon or yield. Where the same data are used in both cases, the correla- 
tion between yield and spectral data appears to be as high as that between 
crop condition and the spectral data (4-11 to 13, 5-28). This similarity 
suggests that the spectral data may contain more yield "information" than 
the Individual crop condition parameters. Analytical definitions or statis- 
tical proofs of the crop condition parameters which are actually "viewed" by 
Landsat data apparently are not available, and more than one crop feature may 
be associated with the yield effects. The yield-Landsat Indicator correla- 
tions also show changes with crop calendar (3-3, 5-27) with the peak correla- 
tion apparently near heading. Thus seasonal, geographical, and culturally 
Induced variability in crop calendars must somehow be taken Into account in 
yield models using Landsat. 

The concept of an advanced yield model consisting of both spectral and 
meteorological components was endorsed (2-17, 4-25, 5-22, 8-4). Rationale 
for using meteorological parameters originates from known between season and 
near harvest dynamics in crop environmental -condition-yield relationships. 

On the other hand, MSS spectral components could both simplify and make 
advanced yield models more accurate by accounting *or the multitude of yield 
affecting factors Integrated by the crop up to the reflectance observation 
time. Studies with the Infrared bands planned for Landsat C and Landsat D 
indicate that they will provide even more crop condition and yield infor- 
mation than the current Landsat data (3-4). 
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As Interim or alternatives to the direct use of Landsat data In yield models, 
the briefings presented two Indirect approaches where Landsat data Is used In 
conjunction with available yield models. Procedures Indicating some feasi- 
bility were those where Landsat estimated LAI Is Input to ET or Growth Models 
(6-13), Landsat Interpretated overrides to crop diagnostic submodels (7-2), 
and monitoring the areal extend of drought (5-1 to 9). The Improvement rela- 
tive to the equivalents of these procedures without Landsat data has not been 
experimentally (statistically) evaluated. 

An Initial list of technical Issues relating to the development of advanced 
yield models was presented (2-18), and several new ones were identified 
during the briefings: 

a Obtaining or Identifying accurate estimates of true yields for specific 
fields (4-14, 5-22, and 8-4). 

e Accounting for within field variability of crop condition, yield, or soil 
moisture for associated Landsat or thermal data (3-6 and 5-22). 

• Separation of crop development stage effects from crop condition variability 
In model development. 

The lack of definition of crop condition or features sensed by MSS data and 
models relating ground observations of these variables to yield could also be 
considered a technical Issue. 

Thus technical rationale and indicators (spot correlations) are the principal 
criteria demonstrating the feasibility of using Landsat data for yield esti- 
mation. To quantitatively assess the feasibility relative to conventional 
methods of estimating yield still requires a considerable amount of model 
development and testing on Independent data. 
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SECTION 2 


LAC IE EXPERIENCE AND OVERVIEW 
OF JSC YIELD PROGRAM 


J. D. Erickson 



LACIE PHASE I AND II ACCURACY SUMMARY 


PRODUCTION 

• WINTER WHEAT 

t U.S, AND USSR ESTIMATES SUPPORTED 90/90 AT-HARVEST CRITERION 

• SPRING WHEAT 

• TENDENCY TO UNDERESTIMATE DUE TO ACREAGE UNDERESTIMATION IN 
U.S. AND CANADA 

• UNDERESTIMATE NOT OBSERVED IN USSR 

YIELD 

• SUPPORTED 90/90 CRITERION IN PHASE I AND II OPERATIONS 

• LOCAL PROBLEMS OBSERVED IN AREAS OF EXTREME WEATHER TESTS 

• 10-YEAR TESTS INDICATED PERFORMANCE MARGINALLY SUPPORTED 90/90 IN USGP 


TECHNICAL MODIFICATIONS IN YIELD FOR PHASE 111 

• INITIALLY, PHASE II TECHNOLOGY WAS USED 
t CCEA I YIELD MODELS ARE IN PHASE II 

• YIELD 

• MODIFICATION OF CCEA I MODELS - IMPLEMENTED APRIL 77 

• EXPANDED TO PREVIOUSLY UNMODELED AREAS IN U.S,, USSR 

• REDEFINED MODEL BOUNDARIES IN U.S. TO ELIMINATE BIASES DUE TO OVERLAP 

• EVALUATION OF SECOND GENERATION YIELD MODELS IN LIMITED AREAS 
i KANSAS, NORTH DAKOTA, 1 USSR WW AND 1 USSR SW OBLAST 

• PRODUCTION 

• INCORPORATED GENERAL ASSESSMENT OF CROP CONDITION BASED ON CLIMATIC 
AND LANDSAT DATA INTO REPORTS 



LACIE PHASE III RESULT5-T0-DATE SUMMARY 


ESTIMATE ACCURACIES 


IN THREE GLOBAL CROP YEARS, LACIE CROP SURVEY TECHNOLOGY HAS PRODUCED 
SIGNIFICANTLY IMPROVED WHEAT PRODUCTION INFORMATION 

• U.S. AND USSR WINTER WHEAT SURVEY ESTIMATES SUPPORTIVE OF 90/90 
CRITERION 1-1/2 - 2 MONTHS PRIOR TO HARVEST 

* 

• PHASE III MODIFICATIONS PRODUCED SIGNIFICANTLY IMPROVED EARLY SEASON 
SPRING WHEAT ESTIMATE IN COMPARISON TO PHASE II - HOWEVER, KEY 
TECHNICAL ISSUES REMAIN WITH SMALL FIELDS/REGISTRATION 

• YIELD ESTIMATES SUPPORTIVE OF 90/90 - TEST AND EVALUATION OF MODELS 
MODELS PLUS POOR PERFORMANCE IN OTHER-THAN-NORMAL WEATHER CONDITIONS 
INDICATES NEED FOR FURTHER IMPROVEMENT 



NASA/JSC YIELD RSD OBJECTIVES 


GENERAL - DEVELOP IMPROVED TECHNOLOGY TO PREDICT MORE ACCURATELY AND WITH KNOWN 
CERTAINTY, VALUE OF YIELD PER HARVESTED ACRE WHICH CAN BE USED IN 
PRODUCTION FORECASTING FOR LARGE U.S. AND FOREIGN REGIONS AT REGULAR 
INTERVALS PRIOR TO HARVEST 

SPECIFIC - OBTAIN MODELS WHICH: 

• ARE UNIVERSAL IN APPLICABILITY WITH A MINIMUM OF ANCILLARY DATA 

» ARE MORE RESPONSIVE TO WEATHER -- ESPECIALLY ABNORMAL AND 
EPISODIC WEATHER 

• INCORPORATE DIRECT OBSERVATION OF CROPS, WEATHER, SOILS, AND SOIL 
MOISTURE FROM SATELLITES 

• SUPPORT THE IDENTIFICATION OF PARTICULAR CROP IN REMOTE SENSING 
PROCEDURES 

• FLEXIBLE AND EFFICIENT TO OPERATE AND UPGRADE 



FY78 YIELD-RELATED SUPPORT' *G RESEARCH & TECHNOLOGY 


t KSU - FEYERHERM 

CONTINUE DEVELOPMENT AND TESTING OF KSU YIELD MODEL 

• KSU - KANEMASU 

CONTINUE DEVELOPMENT AND TESTING OF KANEMASU GROWTH AND YIELD MODELS 

• USDA/ARS - WHEAT YIELD MODELING TEAM 

COLLECT DATA AND DEVELOP AN ADVANCED WHEAT YIELD MODEL 

• USDA/ARS - BLACK (FY77) 

DEVELOP A MODEL OF WINTERKILL PERCENT 

• DPRA (FY77) 

DEVELOP A WINTER WHEAT STARTER MODEL AND IMPROVED CROP CALENDAR MODE! 

DEVELOP INPROVED TECHNIQUES FOR USING METSAT INFORMATION TO INTERPOLATE 
PRECIPITATION AMONG FIRST ORDER STATIONS AND DERIVE SOLAR RADIATION 
ESTIMATES: PROVIDE IN A FORMAT WHICH NOAA CAN IMMEDIATELY USE 


DEVELOP A HYBRID YIELD MODEL BASED, PERHAPS ON THE FEYERHERM AGROMET 
MODEL AND INCORPORATING LANDSAT - DERIVED (ET/ETP) OVERRIDES 


DEVELOP A MET-BASED MODEL WHICH PREDICTS CROP LANDSAT SPECTRAL 
APPEARANCE ON ANY CALENDAR DATE 


ro 

i / 
VJl 


t RFP 


• RFP 


• RFP 




NASA S 77-12053 


YIELDS ARE DETERMINED BY THE INTEGRATED EFFECTS OF 
BASIC SOILS, CLIMATE AND CULTURAL FACTORS BY CROPS 


YIELD MODELLING OBJECTIVE 


ro 

i 

*^i 



FOR EXAMPLE: 

• CROP CALENDER x DISEASES 


• SOIL SOIL VARIETY 

FERTILITY MOISTURE 



OUTPUT SPACE SCALE, m 


YIELD MODEL TYPES 



TIPE SCALE OF INPUTS 
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TEST AND EVALUATION RESULTS - SUMMARY TO DATE 


• BASIC AGROMET REGRESSION MODELS (CCEA) CONDITIONALLY MEET 90/90 HYPOTHESIS 

• POWER OF 10-YEAR TEST RELATIVELY LOW AND INDEPENDENCE OF SAMPLE RESTRICTED 

• OBSERVED AMPLITUDE ABOUT MEANS OF PREDICTED YIELDS IS SMALL COMPARED TO ACTUALS 

• CASES OF BIAS OBSERVED IN INDIVIDUAL MODEL REGIONS 

• CASES OF LARGE SINGLE YEAR ERRORS IN INDIVIDUAL STRATA ARE OBSERVED 

• VARIANCE COMPUTATIONS VALID 

• POTENTIAL IMPROVEMENT BY SECOND GENERATION MODELS INDICATED BUT NOT YET 
DEMONSTRATED ON FULL REGION (90/90) CRITERIA 
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KANEMASU (KSU) 
YIELD MODEL CONCEPT 




ADVANCED SPECTRAL 
AND 

AGRO-METEOROLOGICAL 
YIELD MODEL CONCEPT 
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RATIONALE FOR METEOROLOGICAL-SPECTRAL HYBRID MODELS 
• FOR METEOROLOGICAL PARTS 

• ENVIRONMENTAL CONDITIONS AT CRITICAL TIMES IN CROP'S LIFE CYCLE CAN AFFECT 
YIELD WITHOUT CHANGING ITS APPEARANCE 

• IN MOST GRAIN CROPS THE YIELD COMPONENT IS HIDDEN FROM VIEW AND VARIES IN 
PROPORTION TO THE VISIBLE VEGETATIVE COMPONENTS 

• CAN UTILIZE AVAILABLE DATA SOURCES WHICH ARE INDEPENDENT OF CLOUD COVER 

• INCORPORATE EXTENSIVE MODELING EXPERIENCE 

t FOR SPECTRAL PARTS 

• STANDING CROPS INTEGRATE MOST ENVIRONMENTAL AND CULTURAL EFFECTS UP TO ANY 
POINT IN TIME, THUS HAVE POTENTIAL FOR DIRECT ASSESSMENT OF Nr STRESS OR 
DAMAGE... YIELD 

• CORRELATION OF SPECTRAL DIFFERENCES WITH FIELD-TO-FIELD YIELD DIFFERENCES 
DEMONSTRATED 

• DATA IS NOT SUBJECT TO ACCURACY DIFFERENCES BETWEEN REGIONS 



HNICAL ISSUES 


• DEFINITION OF OPTIHUN AGROPHYSICAL STRATA FOR YIELD AND AREA ESTIMATION 

• QUALITY CONTROL AND STANDARDIZATION OF MODEL BUILDING, TESTING, AND OPERATIONAL DATA 

» SAMPLING ENVIRONMENTAL FACTORS FOR FULL STRATA PREDICTION VS. PREDICTING YIELD FOR A 
SAMPLE OF A STRATUM 

• OPTIMUM COORDINATES AND FREQUENCY FOR INTERPRETATION OF METEOROLOGICAL INPUTS FOR ANY 
PARTICULAR MODEL 

• DEFINITION OF CROP FEATURES VIEWED BY REMOTE SENSORS AT VARIOUS TIMES OF SEASON 

• ACCOUNTING FOR WITHIN AND BETWEEN STRATA VARIABILITY OF CROP CALENDARS, SOIL CHARACTERIS- 
TICS AND MANAGEMENT PRACTICES 

• ADEQUATE AGROMET AND SPECTRAL METHODS TO PREDICT PUNTING AND CROP DEVELOPMENT STAGES 
» ADEQUATE MODELS TO TRACK SOIL MOISTURE AVAILABILITY TO CROPS 

• DEVELOPMENT OF INTERCHANGEABLE SPECTRAL OR METEOROLOGICAL INPUTS FOR A COMMON YIELD MODEL 

• ESTIMATION OF FACTORS WHICH CAUSE TRENDS IN YIELD 

• DEFINING EFFECTS OF EPISODAL EVENTS WHICH ARE NOT TAKEN INTO COMPLETE ACCOUNT BY A MODEL 

• USE OF LANBSAT-C THERMAL BAND DATA IN PREDICTING YIELD 

• UTILIZATION OF SERIAL CORRELATIONS (TIME, SPACE, AND CROP TYPES) INTO YIELD PREDICTIONS 

• APPROPRIATE TEST AND EVALUATION 



SECTION 3 


ON WINTER WHEAT YIELD FROM LANDSAT AND LANDSAT 
FOLLOW-ON SATELLITES 


J. C. Harlan, Jr. 
Remote Sensing Center 
Texas A$M University 
College Station, Texas 


The feasibility for determining winter wheat yield from earth ob- 
servation satellite data has been examined. The desirability of 
utilizing crop observations for yield estimation is due in part to 
the limitations of meteorological yield models. Models based on 
meteorological data are adversely affected by the sparseness of 
weather stations. In addition, yield reducing factors which are 
not weather-related -- such as nsects, diseases and soil fer- 
tility -- are difficult to quantify for inclusion in the models. 


The hypothesis that wheat yield can be determined from multitemporal 
Landsat data is based on work relating grain yield to the size and 
duration of the crop photo synthetic system and on studies relating 
Landsat data to green biomass. Cereal crop grain yield can largely 
be attributed to the photosynthesis during the growth and maturation 
of the grain [1] . The amount of photosynthesis depends on two 
factors: the size and duration of the photosynthetic system; and 

the efficiency of that system. 


The correlation between green biomass and the Transformed Vegetation 
Index has been established [2, 3, 4], The size of the photosynthetic 
system of wheat is reflected in the Landsat measurement of green 
biomass. Likewise, the duration of the system is determinable from 
repetitive Landsat coverage. The system efficiency is variety de- 
pendent, not measurable from Landsat, and, therefore, a noise factor. 


In this study data has been analyzed for selected locations in the 
southern Great Plains region of the United States, from four crop 
years. High resolution spectral data acquired of commercial wheat 
fields were used to simulate data from Landsat and projected Landsat 
Follow-On sensors. Actual Landsat data were used as well. Field- 
by-field yield data were acquired from farmers in terms of actual 
harvested grain weight or from the U.S. Department of Agriculture 
(USDA) in terms of yield/area extrapolated from field samples. 

All Landsat data utilized in this study were treated under a stand- 
ard procedure. The data preprocessing consists of 1) application 
of a cosine correction for sur. angle so that the sun appears to 
have been at zenith; 2) grouping pixels by individual land units 
(in this case, farmers’ fields); and 3) calculation of the Landsat 
band mean vector and correlation matrix for each of the fields. 

The Transformed Vegetation Index (T r !) is calculated from the 
following equation: 


Work supported by the U.S. National Aeronautics nd Space Adminis- 
tration through contract NAS9-14470. 
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TVI - 


Vi 


MSS 7 

itssr 


- MSS5 

♦ MSS5 


+ 0.5 


Where the MSS values are the mean radiances (sun angle corrected) 
for the given field. 


S-191H spectral reflectance data (field spectra referenced to a 
barium sulfate coated panel) were processed to the bandpasses of 
the Landsat MSS plus two other near infrared bands, one of which 
is equivalent to the 1.55-1.75 pm band proposed for the Thematic 
Mapper for Landsat D. The_ bands are given below. Vegetation paran 
eters (VP) of the same form as the ^VI were utilized with the' S- 
191H reflectance data. These parameters each use two of the band- 
passes at one time as follows: 


VP2i 


VI 


and i - Band 2 
Band i + Band 2 


♦ 0.5 


ORIGINAL PAGE IS 

or p- ■ ' rv 


Where i * 3, 4, 5 and 6; the Band values are the mean reflectance 
values for the given field and the band limits are: 


1 0.50-0.60 pm 

2 0.60-0.70 pm 

3 0.70-0.80 pm 


4 0.80-1.10 pm 

5 1.15-1.30 pm 

6 1.45-1. 7S pm 


Preliminary to determining the relationship of yield to spectral 
data the correlation between Landsat TVI, calculated from sun angle 
corrected radiance values, and the S-191H vegetation parameters, 
calculated from reflectance was determined. Landsat-1 and the S- 
191H acquired data on the same date twice during the spring of 1974 
(3/16 and 5/27) and once within a day of each other under the same 
atmospheric conditions (4/3 and 4/4). To compare the data sets from 
the two sensors, the TVI and VP2i values were calculated for all 
wheat fields observed by both. Regression analysis of the TVI ver- 
sus each of the S-191H vegetation parameters produced the results 
expressed in Table 1. Correlation coefficients for TVI versus VP2i 
are given for each date. For each date VP23, VP24_and VP25 are 
seen to be very highly correlated to the Landsat TVI. VP26, where 
band 6 is the 1.55-1.75 pm Thematic Mapper band, is radically differ- 
ent, however, particularly during ripening (5/27). VP26 does not 
measure the same thing, then, and the 1.55-1.75 pm band pass must 
contain unique information. 

The first efforts in the investigation were applied to data collec- 
ted at one site during the 1973-1974 crop year. Seven Landsat ac- 
quisitions of a 4.8 by 4.8 km (3 by 3 mile) commercial farming area 
were examined which covered all growth stages from fall establish- 
ment to ripe. TVI values were determined for each Landsat pass for 
23 wheat fields for which yield data were available. The yield 
values ranged from 0.86 to 4.035 metric tons/hectare (12.8 to 60.0 
bushels/acre). Linear regressions were applied with yield as the 
dependent variable and TVI or (TVI) 2 for each date as the indepen- 
dent variables. The best combinations of variables from thiee or 
less Landsat acquisitions are given in Table 2. More variables did 
not significantly improve the relationship. Results for four or 
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TABLE 1 Correlation Coefficient 
Between S-191H Vegetation Parameters and LAND SAT TVI 


Date 

VP23 

VP 2 4 

VP 2 5 

VP 2 6 

3/16/74 

Jointing 

0.944816 

0.946827 

0.944610 

0.562205 

4/4/74 

Jointing 

0.903736 

0.923658 

0.910313 

0.656690 

5/27/74 

Ripening 

0.899234 

0.882267 

0.854460 

-0.027195 


less passes are the most useful since it is unlikely that more bio- 
phases than that could be consistently acquired of a given site by 
Landsat due to cloud cover or snow on the ground. 


TABLE 2 Bushland 1973/1974 Yield Estimation 
Regression Models from Landsat TVI Values 


R 2 

Number of 
Variables 

Number/Names 
of Biophases 

Order of 
Equation 

0.874 

3 

2/Fall establishment (11/28/73) 
Heading (5/8/74) 

2 

0.861 

3 

2/Tillering (12/16/73) 
Heading (5/8/74) 

2 

0.888 

5 

3/Fall establishment (11/28/73) 
Tillering (12/16/73) 

Heading (5/8/74) 

2 


The results from the Bushland analysis were very encouraging. The 
work since that time has been dedicated to testing the Bushland 
technique. The first data set used in the testing is that from 
S-191H observations over western Kansas farms during 1974/1975. 

From the correlation between TVI and VP24 the ability to simulate 
TVI data from S-191H values was established. For each of the seven 
dates, representing five biophases the simulated TVI values were 
calculated and linear regressions again run, this time with TVI- 
cubed and TVI raised to the fourth power as additional variables. 
The actual biophases and variables chosen were somewhat different 
from the Bushland set. The R2 values were also 1 0% - 20% lower for 
western Kansas. Even so, the results indicate that the Landsat - 
derived parameters by themselves could explain most of the varia- 
tion in the observed yield data. 

The other S-191H vegetation parameters were also calculated for 
each wheat field. Regression analyses were run for each vegetation 
parameter using the VP value and its square as variables to deter- 
mine the capability for explaining yield variation through other 
vegetation parameters besides TVI (VP24) and, specifically, to 
examine the use of the Thematic Mapper band, 1.55-1.75 ym (S-191H 
band 6) . 

Table 3 summarizes the results of the regression analyses utilizing 
the S-191H vegetation parameters. Using no more than four passes 
it is apparent that each vegetation parameter is capable of explain 
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ing most of the yield variation. As noted above the VP26 contains 
different information than VP24 and VP25. This comes through again 
in that different biophases were chosen when using VP26, while the 
resulting R2 values were comparable to those obtained using the 
other two parameters. It can be concluded that although no apparent 
increase in yield estimation accuracy occurred by using the 1.S5- 
1.75 ym information in this vegetation parameter the flexibility of 
using the satellite-borne Thematic Mapper in yield estimation will 
be greater since more combinations of three or four cloud free 
passes could be used with the same accuracy. For example, if til- 
lering and ripening were the only cloud free passes, then TVI 
could be used; while if jointing and heading were good, but tillering 
not, then *'TVI26" could be used where the 2 and 6 were, respectively, 
the MSS 5 equivalent and the 1.55-1.75 ym band. 


TABLE 3 Regression Model Results From Landsat Band 
and Proposed Sensor Band Parameters 


Parameter 

R 2 

Number of 
Variables 

Number/Names 
of Biophases 

Order of 
Equation 

VP24 

0.756 

4 

3/Tillering 

(3/20/75) 

Ripening 

(6/2/75) 

Ripening 

(6/9/75) 

2 

VP 2 5 

0.743 

4 

4/Tillering 

(3/20/75) 

Ripening 

(6/2/75) 

Ripening 

(6/9/75) 

Ripening 

(6/17/75) 

2 

VP26 

0.654 

4 

3/Jointing 
(4/8/75) 
Green Headed 
(5/21/75) 
Ripening 
(6/17/75) 

2 


0.738 

5 

4/Jointing 
Green Headed 
Ripening (6/9) 
Ripening (6/17) 

2 


Results of work relating temporal series of vegetation parameter values 
to wheat yield pointed out that some important yield influencing 
factors were not accounted for in the data set. This conclusion was 
evident because algorithms developed for one location or one year were 
not accurate when applied to another location or even a second year 
at the same location. Consequently two approaches for further work 
were taken. 
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One to be done in-house at NASA/Johnson Space Center involved 
examination of the TVI yield estimation technique in a year-to- 
year mode. That is, determining if estimates of year-to-year 
variation of yields corresponded to variation of TVI at certain 
growth stages. The approach involved combining the temporal series 
of TVI values per year per field or CRD (crop reporting district) 
unit into terms ’TVI (75-76)’ and 'TVI (76-77) ' . These would be 
combined with the known yield for crop year '75-^76, ’¥(75-76)', 
to determine an estimate of the ’76-*77 yield, *Y(76-77)*: 

*(76-77) - YC7S-76) . £$ ^7 # 

The accuracy of these yield estimates would be determined by com- 
paring against ASCS- and SRS-reported yields for '76-’77. 

The second approach, implemented at Texas A§M University, involves 
examining the response of multispectral scanners to occurrences of 
yield detractant phenomena such as drought and disease. The empha- 
sis is being placed on quantifying the relationship between crop 
condition and scanner parameters. The ultimate goal is to increase 
the universality cf Landsat -based yield estimation techniques. 

Johnson Space Center Agricultural Field Measurements Program data 
are used exclusively since there are no other sources of such infor- 
mation. NASA helicopter-borne spectrometer (S-191H) and truck- 
mounted spectrometer data have been used to simulate Landsat MSS and 
Thematic Mapper band values. These are used in conjunction with 
agronomic ground data acquired by the USDA/ASCS and NASA/JSC in 
support of the Field Measurements Program flights and Landsat passes. 
Data have been used from the 1974-1975 and 1975-1976 crop years at 
both the Williams County, North Dakota, Intensive Test Site and 
Agricultural Experiment Station and from the Finney County, Kansas, 
Intensive Test Site and Agricultural Experiment Station. 

At the writing of this document analysis of the responses of indi- 

vidual scanner bands to yield/growth detractant occurrences is nearing 
completion. Analysis of the visible and reflective infrared band 
values substantiates the unique character of the proposed 1.55- 
1.75 urn Thematic Mapper band 5 mentioned above. When fields were 
grouped by ASCS-reported detractant (no detractant, drought, uneven 
stand, or weeds) for a given growth stage the reflectance within the 
Thematic Mapper band 5 was generally different for "detractant" 
groups than for "no-detractant" groups (control fields) . This was 
also true, but to a lesser extent, with a band between 2.10 and 
2.35 um. It was not the case with either the present Landsat MSS 
bands of the Thematic Mapper bands 1 through 4. The Thematic Mapper, 
therefore, appears to afford future analysts the opportunity to 
classify wheat into condition classes or groups, related to probable 
yield, on each satellite pass. 

A separate analysis scenario of the simulated Landsat C and Thematic 
Mapper thermal infrared band data for the Kansas Intensive Test Site 
has been applied. After eliminating advection from upwind fields as 
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a source of radiative temperature difference, within-field causes 
for differences were hunted. Spots within several fields were noted 
as relatively warmer or cooler than the mean field value on not one, 
but two separate dates three weeks apart (April 18 and May 6, 1976). 
A comparison of soil maps and relative canopy density (from aerial 
photos) has shown that, in general, the warm spots occur on patches 
of silt loam soils supporting less dense vegetation than the major- 
ity of the field area. These patches are more droughty than the 
most prevalent soil in the site. On the other hand the cooler spots 
are found in patches of clay which tend to retain moisture. The 
primary response of the thermal infrared band appears to be from the 
combination of moisture and canopy cover. Analysis is proceeding in 
this area in an attempt to qvantify the thermal band sensitivity to 
these scene factors. The thermal band of Landsat C, to be available 
in 1978, should be valuable in interpreting crop condition from 
scanner data. 


1. P. J. Welbank et al.. Annals of Botany, N.S. 30 (118) 291 (1966). 

2. J. W. Rouse, Jr. et al., in: Third ERTS-1 Symposium , NASA- 
Goddard Space Flight Center, 1973, p. 309. 

3. R. H. Haas et al., in: Earth Resources Survey Symposium , NASA- 

Johnson Space Center, 1975, p. 43. 

4. D. W. Deering et al., in: Tenth International Symposium on 

Remote Sensing of Environment. University of Michigan, 

1975, p. 1169. 
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• CRTS— 1 • 

• # 

•SITE PROCESSING REPORT • 

6 ***« «**•*#**•* *#***»•**»* 

SITE DESIGNATOR: GPSO*IT CELLS 2309 TO 2*28 

IMAGE IDENTIFIERS 1619-16530 LINES 200 TO 280 

OATE ACQUIRED: 03APR74 TOTAL POINTS 40 

•••CORRECTED FOR SUN ELEVATION 49 DEGREES”* 


RAOtANCE I MUATTS/SQCM-STR-MI CRON ETERI 



MEAN 

STANDARD 

UAVELENGTH 



DEVIATION 

(MICROMETERS! 

BAND 4 

7.58 

0.54 

.5 - .6 


BANO 5 

5.70 

0.80 

•6 • • 7 


BAND 6 

9.10 

0.98 

•7 • *8 


BANO 7 

9.36 

1.41 

• 8 -ia 



NORMALIZED COVARIANCES 



BANO 4 

BANO 5 

BAND 6 

BANO 7 

BAND 4 

1.000 

0.808 

-0.094 

-0.163 

BANO $ 

0.808 

1.000 

-0.395 

-0.484 

BANO 6 

-0.094 

* -0.395 

1.000 

0.937 

BANO 7 

-0.163 

-0.484 

0.937 

1.000 


BAND RATIO PARAMETER 0.243 fferf •/ TVl) 
TRANSFORMED PARAMETER 0.862 ( TYl, ) 


16.0 ♦ 


8.0 ♦ 


* 





• 0.5 0.6 0.7 0.8 0.9 1 . 0 * 

• RAO! ANCE VS HAVEL DIGTH • 

• PREPAREO BY: TEXAS ACM UNIVERSITY, REMOTE SENSING CENTER • 

• OATA ANALYSIS LABORATORY 22JAN75 • 

• 6 
**«*****•♦*•« *****••*«•*•** ft******************************** 

Figure 2-1. A Typical LANDSAT MSS Data Site Processing Report. 
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LANDSAT AND OTHER PROPOSED SENSOR BANDS 


LANDSAT plus 2 others 

1 0.50-0.60 

2 0.60-0.70 

3 0.70-0.80 
A 0.80-1.10 

5 1.15-1.30 

6 1.45-1.75 


THEMATIC MAPPER 

0.48-0,53 

0.53-0.60 

0.62-0.68 

0.74-0.91 

1.55-1.75 

THERMAL IR 



I 


FINNEY CO. 

REGRESSION MODELS FROM SIMULATED TVI VALUES 


R2 

[ 

NUMBER OF 
VARIABLES 

number/names 

OF BIOPHASES 

ORDER OF 
EQUATION 

: 0.767 

4 

3/til’uERing (3/20/75) 
RIPENING (6/2/75) 
RIPENING (6/9/75) 

3 

0.703 

4 

2/tillering (3/20/75) 
RIPENING (6/9/75) 

4 

0.777 

4 

3/tillering (3/20/75) 
ripening (6/2/75) 
ripening (6/9/75) 

r 

4 

• — _ 
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CONCLUSIONS FROM PREVIOUS CONTRACT WORK: 

1) High correlations were achieved between wheat yield and 

MULT I TEMPORAL VALUES OF LANDSAT MSS AND ThEMATIC MaPPER 

vegetation indices for each individual location and/or 

YEAR. 

2) Regression models differed/ however, for different locations 
and/or years. Factors not. accounted for by the vegetation 
parameters need to be examined. 

3) Conclusion 2) led to the present contract work as well 

AS A JSC in-house project to estimate year-to-year variation 
in yield from corresponding variation of the vegetation 
parameter TVI at selected growth stages. 

4) a non-destructive LAI (leaf area index) measurement technique 
was developed from photography and LAI MEASUREMENTS taken 

for the Agricultural Field Measurements Project at the 
Finney County supersite. 
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TAMU/RSC 


ISSUE ; IMPROVED YIELD MODELS 

PRESENT TASK ; STRESS, STAND QUALITY AND CROP CONDITION FROM 
SCANNER DATA 

OBJECTIVES RELATED TO ISSUE ; 

1) Determine the response of multi spectral scanners 
(Landsat 1 and 2 MSS, Landsat C MSS and Thematic 
Mapper) to occurrences of yield detractants such 
as drought and disease. 

2) Develop algorithms relating spectral response to 
yield. 

APPROACH : 

1) Use FSS and truck spectrometer data from Finney County and 
Williams County Supersites for 1974/1975 and 1975/1976 to 
determine Landsat MSS and Thematic Mapper reflectance values. 

2) Calculate vegetation parameters from the reflectance values 

AND CORRELATE THEM WITH CROP CONDITION TO DETERMINE STRESS 
EFFECTS . 

3) Summarize results and analyze and interpret them in terms 
of identification and quantification of stresses and their 
effects. 
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1975 & 1976 KANSAS WINTER WHEAT: 
DEPENDENCIES OF SPECTRAL DATA ON GROUND DATA 


SCANNER BAND WITHOUT DETRACTANTS WITH DETRACTANTS 


HSS 6 

Stand Quality 

Detractants ^ 


WEEDS, 

uneven 

MSS 7 

Stand Quality 

Detractants 


stands 

Thermal 

Canopy Height 

Canopy Height 


(8.0 - 13.5 fin) 





TM 3 

— 

Detractants > 

j 


WEEDS, 

uneven 

TM H 

Stand Quality, 
Canopy Height 

Detractants J 


stands. 

Tfl 5 

Canopy Height 

Detractants, \ 

Canopy Height ( 

1 

drought 

both 

2.10 - 2.35 tin 

Stand Quality 

Detractants 

) 

years 

TVI7 

Canopy Height 

Detractants, 





Canopy Height 



•Growth stage and percent ground cover were both correlated with 

EACH BAND. 


Stand Quality - 1 to 6 

Growth/Yield Detractants - coded 

Growth Stages - 1 to 10 

Ground Cover - 1 to 5 in 20% increments 

Canopy Height - in inches 
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Percent Reflectance 


Winter Wheat, Kansas 1975. Growth Stage: Beginning to head 



MSS4 MSSS MSS6 MSS7 TM5. 2 n 2pm 


Growth/Yield 

Detractant. 

Q No detractant 

* Drought 

* Uneven stand. 


Multispectral Scanner Bands 

YIELD/GROWTH DETRACTANT CLASS MEAN REFLECTANCE VALUES IN SELECTED SCANNER BANDS 



ercent Reflectance 


! Winter Nheat, Kansas 1975. Growth Stage: Beginning to head 



YIELD/GROWTH DETRACTANT CL ASS MEAN REFLECTANCE VALUES IN SE 


' Growth/Yield- 
( Detractants 

No detractant* 

* Drought 

i 

* Uneven stand- 


.SCANNER BANDS 



rccnt Reflectance 


Winter Wheat, Kansas 1975. Growth Stage: Fully headed 




Percent Reflectance 


Winter Wheat, Kansas 1975, Growth Stage: Fully headed 



YIELD/GROWTH DETRACTANT CLASS MEAN REFLECTANCE VALUES IN SELECTED SCANNTER BANDS 



Percent Reflectance 


Winter Wheat, Kansas 1976. Growth Stage: Beginning to ripen 




Percent Reflectance 


Winter Wheat, Kansas 1976. 


Growth Stage: Beginning to ripen 


35 
30 - 


25 - 

I 

t 

i 

! 20 - 

I 

15- 
10 - 

5 - 



I 





I 


! Growth/Yield 
| Detractants 
I No detrar.tant 

1 4 ; 

j+ Drought 

I* Uneven stand 
> Weeds 


'"i " ■ " 

TM1 


— l~ i 1 — 

TM2 TM3 TM4 

. Multispectral Scanner Bands 


YIELD/GROWTH DETRACTANTS CLASS MEAN REFLECTANCE VALUES IN SELECTED SCANNER BANDS 


Percent Reflectance 


^Winter Wheat, Kansas 1976. Growth Stage: Ripe, mature 



YIELD/GROWTH DETRACTANT CLASS MEAN REFLECTANCE VALUES IN CERTAIN SCANNER BANDS/ 



Percent Reflectance 


iWin.ter .Wheat, Kansas 1976. Growth Stage: Ripe, mature 



\ 

Multi spectral Scanner Bands 


YIELD/CROWTH DETRACTANT CLASS MEAN REFLECTANCE VALUES IN. SELECTED SCANNER BANDS/ 


ADVECTIVE INFLUENCES ON 
MEAN FIELD TEMPERATURE 


April 18, 1976 


Flight Line 
Number 

Dissimilar Fields 
Upwind . 

Similar Fields 
Upwind 

•10 

82.98* 4 

35.013 1 

9 

83.423 4 

83.342 5 

8 

33.586 5 

82.393 1 


X = 83.416 

X = 83.332 


May 6, 1976 


Flight Line 
Humber 

Dissimilar Fields 
Upwind 

Similar Fields 
Upwind 

10 

73.312 3 

75.427 1 

9 

74.298 4 

73.425 5 

8 

73.913 5 

73.480 1 


X = 73.891 

X = 73.755 


O 1 1 1 

* Units are mw cm l ' sr x 10 . Spectral bandwidth 
is 8.0 - 13.5/<m. 
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COMMON HOT/COLD AREAS 
ON 4/18 & 5/6 


ORIGINAL PACT; r.j 
OF POOH Ou-iLli* 


FINNEY COUNTY; KANSAS (SUPER SITE) 7 
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TARU/RSC 


PUNS FOR FURTHER ANALYSIS : 

• Completion of present tasks by November 30 , 1977 
(end of contract). 

• For "Improved Yield Models" - Develop algorithms 
relating spectral response to yield. 

• For "Lanbsat C Use" - Complete examination of thermal 

BAND SENSITIVITY TO CHANGES IN SCENE FACTORS. DELINEATE 
T..ESE RELATIONSHIPS FOR USE AS INTERPRETATION TOOLS, 

• For "R & P Data Set Availae t lity" - Complete development 

AND TEST THE NON-DESTRUCTIVE UI TECHNIQUE. 
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WORLDWIDE WHEAT PRODUCTION FORECASTS USING LANDSAT DATA* 
Richard F. Nalepka* 


ABSTRACT 

Discussion is presented of tb philosophy, background, and activities 
carried out at ERIK to utilize Landsat data to help forecast the yield and 
production of wheat. Results are presented which demonstrate the empirical 
relationships between wheat yield and percent green wheat cover, percent 
green wheat cover and a Landsat green measure, and wheat yield and the Land- 
sat green measure. Correlations of early season Landsat estimates of yield 
with fanners harvested yield are shown to be as good or better than more 
conventional estimates made later in the growing season. The variance in 
yield accounted for by Landsat variables is also shown to parallel that 
accounted for by several important cultural variables (detailed information 
on these variables would normally not be available in an operational system) 
Results of yield prediction extension are also presented. 

A discussion of a new direct production forecasting procedure using 
Landsat data is presented which potentially overcomes many of the serious 
prcblems (e.g., small fields and cloud cover over specific sites) being 
faced by other available approaches. Initial test results are presented 
which demonstrate quite accurate early season forecasts of production over 
regions as small as LACIE sites and as large as a crop reporting district. 

Further activities are recommended to investigate the use of Landsat 
data for identifying crop condition and estimating yield and to investigate 
the joint use of Landsat data, Metsat data, and agromet models. A strong 
recommendation is made that direct wheat production fnrpracrlng procedures 
should be further developed and evaluated. 

- 

Presented the Landsat Crop Condition and Yield Briefing held at 
NASA Headquarter^ on September 27, 1977. 

** 

Mr. Naleptca is the Head of the Multispectral Analysis Section of the 
Environmental Research Institute of Michigan's (ERIM) Infrared and Optics 
Division . . 
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WORLDWIDE WHEAT PRODUCTION FORECASTS USING LANDSAT DATA 


Principal Investigator 

Richard F, Nalepka 

Co- I nvest i gator 

John E. Colwell 


Presented at the Landsat Crop Condition and Yield Briefing 

NASA Headquarters 
September 27 , 1977 



OUTLINE 


• Basic Philosophy 

• Background 

• Description of Sites Examined 

• Types of Activities 

• Issues Addressed 
•Discussion of Results 

• Conclusions 

• Recommendations 



BASIC PHILOSOPHY 

i 

i 

I 

! At any point in time the crop itself best represents and integrate 

I THE EFFECTS OF VARIABLES SUCH AS: 

j 

! - Planting Date 

; - Available Sunlight 

! - Available and Useful Moisture 

i - Hail or Wind Damage 

- Winterkill 

- Fertilization i 

' - Insect and Disease Damage 1 

I 

\ - Farming Practices 




BACKGROUND 


• Theoretical study for NASA/JSC in early 70' s using ERIM growth and 

' CANOPY REFLECTANCE MODELS TO INVESTIGATE THE POSSIBILITY OF 
SUCCESSFULLY USING SATELLITE MSS DATA TO AID IN FORECASTING 
WHEAT YIELD 

I 

I 

i • Empirical investigation for NASA/GSFC to establish how well? 

- Wheat yield is related to field vegetative condition 

- Landsat data can be used to estimate field vegetative 

condition 

- Landsat data can be used to help forecast wheat yield 

(production) 




SITES EXAMINED 


KANSAS (74-75 & 75-76) 

• LACIE Intensive Test Sites 

- Finney (old and new) 

- Ellis 

- Rice 

- Sali.me 

i 

• LACIE Blind Sites in Central Crop Reporting District 

• Central Crop Reporting District 


£Fmm 



LOCATION OF WHEAT FORECAST TEST AREAS 
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ACTIVITIES 

• Field Measurements and Sample Collection 

• Laboratory Measurements 
•Data Reduction and Analysis 
•Model Calculations and Analysis 

• Lands at Data Processing and Analysis 



ISSUES ADDRESSED IN SATISFYING INVESTIGATION OBJECTIVES 


• Optimum Single Time (near heading) 

• Landsat Green Indicators 

• Comparison With Alternate Approaches 

• Importance of Various Cultural Factors 

• Model Extension (Geographically and Temporally) 

• Data Screening 

• Direct Production Forecasts 
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LANDSAT GREEN INDICATORS EXAMINED 


• Tasseled Cap Green Channel 
: • MSS7/MSS5 = R75 

• VMSS77MSS5 - SQ75 

• V(M$S7 - M$S5)/(MSS7 + MSS5) + 0.5 ■ TVI 

• MSS') - MSS7 + 96 = G 




Wheat Yield 


75.000 + 
+ 

60.000 + 
+ 

45.000 + 
+ 

30.000 + 
+ 

15.000 ♦ 
+ 


* * 


* 



* 


* 


* 


* 


0 . 


+ 

+ + + + + — — + + + + + — — + 

0 . 30.000 60.000 

15.000 45.000 75.000 

Percent Green Wheat Cover 


ERIM MEASUREMENTS OF PERCENT COVER VS WHEAT YIELD 

{Combined 1976 and 1975 Data) 

FINNEY ITS 

2 YRS, 9 FIELDS EACH, HRIM 

MID-MAY ' i ^ 
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75.000 + 



* 


* 

* 

* 

* 


~+ + + + + + 

1.0900 

.99000 . 1.1900 


Measure (SQ75) 


MEASURE VS ERIM MEASUREMENTS 
EN WHEAT COVER 




Wheat Yield 


75.000 


60.000 + 


* * 


* 

* * 


A 5.000 + 

* 

+ 

* 

* * 

30.000 + * * * 

*• * 22 * * * * * 

* * * * * 

•f ****** ** 

*2 .* * 

** 2 * 

15.000 + ** 

** * 

+ * 

. * 

0 . + • 

+ + + + + + + + + + + 

.75000 .95000 4.1500 

.05000 1.0500 1.2500 

Landsat Green Measure (S075) 


LANDSAT GREEN MEASURE VS WHEAT YIELD 
CORRELATION = .8 
1 SITE/DATE 
1 PT/FIELD 

LATE APRIL-EARLY MAY 


Term 

Lmrnk 
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CORRELATIONS OF FARMERS YIELD WITH FIELD ESTIMATES 
AND LANDSAT ESTIMATES OF YIELD 


Yield 


Estimator 

Site A 

Site B 

Site C 

Average 

FCIC* 

0.95 1 

0.26 1 

0.74 1 

0.65 

Stand ## 
Quality 

0.97 1 

0.78 1 

0.89 1 

0.71 

Landsat 

(4 Bands) 

0.94^ 

0.80 4 

0.79^ 

0.84 

Landsat 

(TVI) 

0.93^ 

0.79 4 

0.64 3 

0.79 


Dates when estimators were available: 

^Pre-harvest (mid-late June); ^15 April; ^21 May; ^6 May 

"Federal Crop Insurance Corporation objective estimates. 

## 

Agricultural Stabilization and Conservation Service .subjective 
estimates, 
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PERCENT OF VARIANCE IN YIELD ACCOUNTED FOR 
SEPARATELY BY SEVERAL CULTURAL FACTORS 


Percent of 

Cultural F actors Variance 


Planting Date 


0,1 

Wheat Variety 


10,6 

Fallow Previous Year 

(yes/no) 

35.8 

Irrigation (yes/no) 


56.3 

Fertilization (yes/no) 

55.0 

Amount Fertilization 

(lb/acre) 

57 A 
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PERCENT OF VARIANCE IN YIELD ACCOUNTED FOR BY SEVERAL 
COMBINATIONS OF CULTURAL AND LAHDSAT VARIABLES 



Percent Standard 

Variables 

Variance 

Error 

.1-6 (ALt CULTURAL VARS) 

74,9 

6,89 

7-10 (all Landsat VARS) 

87.3 

4.78 

4,5,7,10 (optimum four vars) 90,7 

4,10 

1-10 (all vars) 

93.6 

3.65 

V 'LL Key: 

s 


1 “• VARIETY 

6 ■ AMOUNT FERTILIZER 

2 a IRRIGATION 

7 - SQ75 (May 6) 


3 = FERTILIZATION 

8 = S075 (June 2) 


4 * PLANTING DATE 

9 - SQ75 (June 1?) 


5 a CROPPING 

10 - SQ75 (April 18) 




TWO TESTS OF EXTENSIONS OF LANDSAT WHEAT YIELD PREDICTION 


Landsat RMS Error 1 


From 

To 

Predictor 

Local 

Non-Local 

Bias^ 

21 May 

20 May 

A Bands 

A. AO 

6.70 

-5,00 

Site A 

Site A 

SQ75 3 

5.2A 

5.08 

0,00 



TVI i * 

5.03 

A. 88 

0.02 

18 April 

18 April 

A Bands 

7.A1 

9.10 

-0.23 

Site A 

Site B 

S075 3 

8.12 

10.18 

2.15 



TVI 1 * 

7.98 

9.29 

1.17 


*On field by field basis/ ^Average difference between actual 

IN ^SI-ILLS, AND PREDICTED YIELD/ IN BUSHELS, 


3 VHSS77fiSS5 



-HSS5 ) / (MSS7+HSS5 ) +0 , 5 


2Trih 
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MAJOR PREMISES OF ERIM DIRECT PRODUCTION FORECASTING 


• At a specific time or times in the growth of winter wheat one 

CAN ESTABLISH A STABLE RELATIONSHIP BETWEEN A LANDSAT GREEN 
MEASURE AND THE PRODUCTION OF WHEAT 

• AS A RESULT OF SPECTRAL DIFFERENCES OR SPECTRAl./TEMPORAL DIFFERENCES 
BETWEEN WHEAT AND NON-WHEAT/ NON-WHEAT PIXELS WILL CONTRIBUTE 
MINIMALLY TO THE FORECASTS OF WHEAT PRODUCTION 

• LANDSAT PIXELS CONTAINING BOTH WHEAT AND NON-WHEAT (E.G., BOUNDARY 
PIXELS) WILL PROVIDE AN INTERMEDIATE GREEN MEASURE THEREBY LEADING 
TO FORECASTS OF PRODUCTION FOR SUCH PIXELS WHICH ARE INTERMEDIATE 
AND CORRECT, 


2Frim 
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ERIM DIRECT PRODUCTION FORECAST PROCEDURE 


• Define production-predictive relationship (based on previous years' data) 

• Stratify region to be processed according to crop calendar 

• Select previously established production-predictive relationship appropriate 
to stratum 

| • Automatically screen landsat data to define bad data, clouds, cloud shadows, 
dense haze, and non-wheat categories such as water, trees, and- urban areas 

• For each stratum and each pixel to be processed (perhaps a sample or 
perhaps all non-screened pixels) determine landsat green measure and 
estimate production 

• Determine final stratum production figure by adjusting accumulated stratum 
production to account for screened pixels 


2£R!M 
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POTENTIAL ADVANTAGES OF ERIM DIRECT PRODUCTION FORECAST PROCEDURE 

• Provides an early season estimate 
Eliminates need to locate and identify fields 

• Provides an approach to operating in regions of small 
or irregularly shaped fields 

• Accounts for non-uniformities in fields 

• Addresses reduced total production due to disease, 
drought, etc, 

• May eliminate need for yearly training 

• Eliminates need for identifying specific sites in advance 

• Potentially reduces effect of cloud cover and sampling error 
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INITIAL RESULT FROM ERIM DIRECT WHEAT PRODUCTION FORECAST PROCEDURE 

(Two LACIE Intensive Test Sites) 


Site 

Landsat 

Overpass 

True 

Production 

ERIM 

Production 

Forecast. 

Error 

(Z) 

A 

S May 76 

AO, 600 bu 

A2,700 bu 

5.2 

A 

18 Apr 76 

AO, 600 bu 

A2,800 bu 

5. A 

B 

6 May 76 

27,900 bu 

2A,700 bu 

-11.5 

A+B 

6 May 76 

68,500 bu 

67,A00 bu 

1.6 





FURTHER RESULTS FROM ERIM DIRECT WHEAT PRODUCTION FORECAST PROCEDURE 

(Ten Counties of Kansas Central CRD) 



2 ITS * 3 BUND SUES 


Landsat 

True 

ERIM Production 

Error 

Overpass 

Production 

Forecast 

(percent) 

i 

17 APR 76 

5,38 x 10® Bushels 

5.2A x 10® Bushels 

2.6 



ERIM Production Forecast 
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CONCLUSIONS 


• Landsat data can be effectively used to estimate certain variables which 

ARE REQUIRED IN EXISTING YIELD MODELS (SUCH AS LAI OR PERCENT COVER) 

; • Landsat indicators of yield are as highly correlated with individual field 

YIELD AS ARE ESTIMATES USING TRADITIONAL FIELD SAMPLING METHODS, EVEN WHEN 
USING LANSAT DATA COLLECTED SEVERAL WEEKS BEFORE THE F , ,LD SAMPLES ARE MADE 

' • A CONSIDERABLE AMOUNT OF THE VARIANCE IN INDIVIDUAL FIELD YIELD WHICH IS NOT 
EXPLAINABLE BY METEOROLOGICAL DATA CAN BE ACCOUNTED FOR BY LANDSAT DATA 

• In order for Landsat data to be of maximal use in an operational system, 

IMPROVEMENTS IN THE ABILITY TO REMOVE THE EXTERNAL EFFECTS (PARTICULARLY 
ATMOSPHERIC EFFECTS) ARE REQUIRED 

• The POSSIBILITY OF MAKING DIRECT WHEAT PRODUCTION FORECASTS USING EARLY- 
season Landsat data looks very promising 




RECOMMENDATIONS 


Investigations of the use of landsat data to identify crop condition 

AND ESTIMATE YIELD SHOULD CONTINUE 

Joint use of landsat data, metsat data, and agromet models should 

BE EXAMINED 

Direct wheat production forecasting procedures should be further 

DEVELOPED AND EVALUATED 



SECTION 5 

MONITORING DROUGHT AND YIELD COMPONENTS BY LANDSAT 
0. R. THOMPSON 

Abstract 

In the Large Area Crop Inventory Experiment, Landsat images and 
digital data were used to detect and monitor the drought that occurred 
in the U.S. Great Plains during the 1976 wheat growing season. Landsat 
color infrared images (100 by 100 nautical miles) were used to determine 
and monitor the areal extent. The drought area was rated subjectively 
as to the acreage affected by comparing the 1976 and 1975 Landsat imagery. 
A technique was devised using a vector transformation of Landsat digital 
data to indicate when vegetation is undergoing moisture stress. A rela- 
tion was established between the remote-sensing-based criterion (the 
Grean Index Number) and a ground-based criterion (Crop Moisture Index). 

Landsat was shown to be correlated to plant properties that influence 
yield. Direct correlation of Landsat to yield appears to be feasible 
only at specific growth stages. The use of Landsat for yield estimation 
is difficult because the biological system is dynamic and because of 
atmospheric effects on Landsat. Some problems exist in the different 
methods of acquiring ground truth (yield estimations) and the variations 
that exist among and within fields. However, assessing yield from 
Landsat appears to be feasible. 
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Introduction 


A Large Area Crop Inventory Experiment (LACIE) has been under- 
taken jointly by the U.S. Department of Agriculture (USDA) , the 
National Oceanic and Atmospheric Administration (NOAA) of the 
U.S. Department of Commerce, and the National Aeronautics and 
Space Administration (NASA) to prove an economical application 
of remote sensing from space (3) . The experiment is being con- 
ducted over three consecutive crop seasons in a 3-1/2-year 
timespan and is divided into three corresponding phases (3) . 

Each phase is designed to build on the experience of the previous 
phase or phases. Phase I was conducted during the 1975 crop year 
and concentrated on a system test to identify and estimate the 
wheat acreage within selected major U.S. wheat growing regions 
and to evaluate wheat recognition analyses in other selected 
regions throughout the world. Phases II and III concentrated on 
bringing all elements of a system together in a quasi-operational 
environment to test the technological capability of developing 
area, yield, and production estimates for U.S. test regions and 
other major wheat producing regions of the world. During Phase II 
(crop year 1976), the drought that occurred in the U.S. Great 
Plains wheat growing area w.^s detected and monitored using 
Landsat data (4,5,6). The approach and results of this study 
are presented in this paper. 

Methods 

Two approaches were devised for monitoring drought using 
remote-sensing-based criteria. One approach utilizes color 
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infrared transparencies of Landsat scenes (100 by 1C: autical 

miles) to determine and monitor the areal extent of drought (4,5). 
The other approach utilizes LACIE sample segments (5 by 6 nautical 
miles) and Landsat digital data to indicate automatically when 
an area is undergoing moisture stress (6) . These two methods 
will be referred to, respectively, as Landsat imagery approach 
and Landsat digital approach throughout the paper. 

Landsat Imagery Approach 

The Landsat imagery approach utilized meteorological data to 
initially locate the area where potential drought might occur. 

Once an area was flagged and delin' ated from meteorological data, 
Landsat color composite transparencies, prepared from band 4 
(0.5 to 0.6 micron-), band 5 (0.6 to 0.7 micron), and band 7 (0.8 
to 1.1 microns) , were used to refine the delineation of the 100- 
by 100-nautical-mile area. These color transparencies were 
evaluated by comparison to Landsat imagery of essentially the 
same date in previous years and also to previous b day acquisi- 
tions of the current year. Normal green vegetation on the ground 
is recorded on the Landsat color composites as a bright red 
color. As moisture stress browns the vegetation on the ground, 
Landsat-recorded signatures correspondently decrease in redness. 
Thus, by relating the lack of redness in the signatures where red 
signatures should be present, the areal extent of the drought 
was monitored and delineated by compiling a mosaic of Landsat 
images over the potential drought area. Within the drought area, 
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the effect of the drought upon the wheat crop was evaluated 
subjectively by comparison with the previous year's Landsat data. 
The area was monitored at 9-day intervals until harvest of the 
wheat crop. 

Results of Landsat .imagery Approach 
U.S. Southern Great Plains 

The drought that occurred in the 1975-76 winter wheat crop area 
originated in the summer of 1975 when the soil moisture supply 
was not recharged after the 1974-75 harvest. This acute moisture 
shortage covered a perir ? of over 30 days, between planting and 
emergence of the wheat. During the 1975 Thanksgiving week, a 
major storm system moved through the Great Plains, bringing 
blizzard conditions to most of the O.S. Great Plains. The com- 
bination of these conditir ns caused the winter wheat to go into 
dormancy with very little root system or top growth. These areas 
were monitored fr^m planting using Landsat imagery. At the 
start of spring greenup, it became apparent that portions of the 
U.S. Great Plains winter crop were affected by the extreme dry 
conditions. LACIE monitored the area every 9 days until harvest. 
The drought-affected area in the U.S. southern Great Plains was 
determined from Landsat to be located in the southwestern corner 
of Kansas, in southeast Colorado, and in the Oklahoma and Texas 
Panhandles. The areal extert of the affected area as of April 1, 
1976, is shown ir. figure 1. The drought severity within the 
area was rated subje *. ^ely by comparing the 1975 and IjIS 
Landsat imagery These ratings corresponded to the acreage 
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losses developed from ground-based observations. The Crop 
Moisture Index (CMI) for April 3, 1976, shows that this general 
area was undergoing moisture stress (figure 2) . 

U.S. Northern Great Plains 

The droughts in the ".S. northern Great Plains also originated 
in the summer of 1975 when subsoil moisture was not fully 
recharged. Precipitation was adequate tor winter wheat from 
emergence to spring greenup. Spring wheat had adequate moisture 
for planting, emergence, and early growth; however, lack of sub- 
soil moisture and spring rains caused moisture stress by mid -May. 

LAC IE, using techniques developed from the O.S. southern Great 
Plains drought study, indicated a potential for drought damage 
in the O.S. northern Great Plains by early May. The areal extent 
of the drought was determined from Landsat full-frame color 
infrared transparencies by monitoring the full-frame Landsat 
images from April 18, 1976, until harvest. 

The initial drought-affected area, as determined from full- 
frame images, was located within South Dakota. From April 18 
to June 12, 1976, the area appeared to be deteriorating, but the 
full-frame imagery did not indicate severe effects. The June 11 
to 13 overpays did show the effects of the drought. The area 
delineated at this time continued to expand until the July 8 to 
11 overpass when the drouoht area stabilized (figure 3) . From 
this overpass- the drought area was rated subjectively as having 
been severely or moderately effected. The July 10, 1976, CMI 
shows that this . was under se*«.re moisture stress (figure 4). 
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Lands nt Digital Approach 

The Landsat imagery approach involves the subjective judgment 
of the analyst- interpreter in deciding that a region is or is not 
drought affected. 

A procedure was devised in an attempt to quantify the sub- 
jective judgment of the analyst- interpreter . The data used were 
Landsat multispectral scanner (MSS) values for LACIE sample seg- 
ments throughout South Dakota, which were acquired during the 
1975 and 1976 crop years (figure 5) . This procedure, which uses 
the remote-sensing-based criterion to detect and monitor crop 
moisture deficiencies without analyzing a long record of clima- 
tological data, was evaluated against the CMI, which is developed 
from ground-based meteorological data. 

This procedure, the Green Index Number (GIN) , was developed 
using ideas presented by Kaurh and Thomas (1) The four Landsat 
channels are rotated into the Kauth and Thomas greenness and 
brightness vectors. Each vector is inspected automatically, and 
any vector having values unreasonable for agricultural data is 
discarded. From these data, a green number is computed. The 
green number indicates the density and vigor of vegetation. Once 
the green numbers are computed for each picture element (pixel) 
within the 5- by 6-nautical-mile sample segment, the GIN is com- 
puted. The GIN then is an estimate of the percentage of pixels 
in a Landsat scene having green numbers high enough (>15) to 
indicate full cover of green vegetation. It is computed using 
only Landsat data. 
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The plot of GIN versus time for a normal, predominantly wheat 
sample segment should follow a curve such as a in figure 6. If 
an observed point for a segment fell into the shaded region, the 
segment was classified as drought affected. The bounds for the 
shaded region were defined empirically as shown in figure 6, with 
t defined as the approximate spring emergence date in days. For 
different areas cr years, the shaded area can be moved from side 
to side to cuz^c.'i the greenup curve. The initial point in South 
Dakota was usually near day 110 (t = 110). rhis classification 
was compared to a classification based on the CMI for a Crop 
Reporting District (CRD) , wherein a CRD was classified as drought 
affected if the CMI fell below -0.5 for 2 consecutive weeks. 

Both classifications were restricted to similar time frames. 
Classification was performed only for data between April 1 and 
July 10. 

Results of Landsat Digital Approach 

The data used in the digital approach study consisted of all 
LACIE sample segments in South Dakota which had at least 5 percent 
wheat as measured by the LACIE Classification and Mensuration 
Subsystem [CAMS (2)] in the 1976 growing season. This definition 
yields 17 segments (figure 5) with 34 possible classifications. 

Of the 34, 4 had either insufficient data during the gro ‘ing 
season or data were inaccessible for other reasons. The final 
data set contained 22 segment jears for 13 LACIE segments 
(table 1}. (NOTE: A segment year is defined as an observation 
of 1 segment for 1 year.) The contingency table (table 2), 
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TABLE X. RESULTS OF GIN AND CMI CLASSIFICATIONS 


(From reference 6] 



— 

1*7S 

197S j 

GIN 

CMC 

GIN 

OX 

A 

W 

« 

- 

D 

• 

- 

w 

m 

O 

C 

w 

w 

o 

O 

D 

w 

w 

D 

o 

e 

w 

IV 

0 

o 

F 

w 

w 

m 

D 

C 

w 

IV 

w 

O 

H 

w 

w 

0 

w 

1 

w 


- 

m 

J 

w 

TV 

D 

o 

K 

w 

w 

0 

o 

L 

0 

0 

0 

D 

M 

- 

D 

m 

O 


TABLE 2. CONTINGENCY TABLE OF GIN AND CXI CLASSIFICATION METHODS 

[From reference 6] 
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which applies the two classification methods to the 22 good 
segment years, shows that the classifications based on the CMI 
and GIN are related. It was concluded that the GIN is detecting 
moisture through crop condition. 

One example of the segment classification procedure is shown 
in figure 7. The GIN indicates that 1975 was normal for the 
entire crop season for segment J. In 1976, the GIN indicates 
that by May 24 there was moisture stress in segment J, which 
indicates that the GIN - Jtected vegetation moisture stress at 
the same time as the CMI. 

Conclusions 

Landsat full-frame color transparencies provide a means of 
locating, delineating, and monitoring areal p't'-ont of moisture 
stress over large areas. A technique was de ed using 
Landsat digital data for 5- by 6-nautical-mile sample segments, 
which indicates when agricultural vegetation is undergoing mois- 
ture stress. A relationship bet^en this technique, which 
utilizes remote sensing, and a ground-based criterion (the CMI) 
has been shown. Indications are that Landsat is capable of 
detecting crop moisture deficiencies in areas of the world where 
ground information is not available or reliable. 
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FIGURE 1 . AREAL EXTENT AND EFFECT OF DROUGHT ON APRIL 1, 3.976. 

[From reference 5.J 
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FIGURE 3. DROUGHT-AFFECTED AREA AS DETERMINED FROM LANDSAT FULL- 
FRAME IMAGERY FOR JULY 8 TO 11 AND JULY 11 TO 20, 1976. iFrom 
reference 4 . ] 
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NORTH OAKOTA 



FIGURE 5. SEGMENT LOCATIONS. Map of South Dakota showing 
locations of LACIE 5- by 6-nautical-mile sample segments. 
[From reference 6.] 



FIGURE 6. PLOT. GIN versus time for a normal, predominantly 
wheat segment. [From reference 6.] , 
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for segment J. [From reference 6.] 
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THE GREEN INDEX NUMBER (GIN) PROGRAM PROVIDES AN AUTOMATIC PROCEDURE 
FOR DETECTING AND MONITORING CROP STRESS OVER LARGE AREAS. 


• AN EXAMPLE OF HOW THIS PROGRAM IS USED IN LACIE IS SHOWN IN 
FIGURES 8-10 

• FIGURE 8 SHOWS THE RESULTS OF THE GIN PROGRAM FOR ONE LANDSAT 
PASS OVER THE USSR SPRING WHEAT REGION. MUCH OF THE AREA WAS 
UNDERGOING MOISTURE STRESS 

• FIGURE 9 SHOWS THE NEXT LANDSAT PASS OVER THE AREA AND INDICATES 
STRESS IS STILL OCCURRING OVER THE REGION 

• FIGURE 10 IS A COMPOSITE OF THE TWO LANDSAT PASSES AND SHOWS 
WHERE MOISTURE STRESS OCCURRED DURING JULY 1977. 
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Figure 8.- Moisture conditions over U .S.S.P. spring wheat from the LAC IE Green Index 
Number (GIN) monitoring program (Landsat data acquired June 23, 1977, and July 2 
through July 19, 1977). 




Figure 9.- Moisture conditions over U.S.S.R, spring wheat from the LACIE GIN 
monitoring program (Landsat data acquired July 19 through August 4 , 1977)* 








Figure 10. - Moisture conditions over U.3.S.R. spring wheat frori the UVCIE GIN 

monitoring program for July 1977. 




• LANOSAT CORRELATION TO YIELD MAY BE DIFFICULT BECAUSE: 
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BIOLOGICAL SYSTEM IS DYNAMIC 

• THE CROP CONDITION WHEN LANDSAT PROVIDES BEST CORRELATION TO YIELD 
MAY NOT BE MAINTAINED TO HARVEST 

• TABLE 1 

• LANDSAT ACQUISITIONS MAY MOT BE SYNCHRONIZED WITH CRITICAL GROWTH STAGE 

VARIATIONS EXIST AMONG AND WITHIN FIELDS 

• SPECTRAL 

• FIELD CENTER PIXELS 
t BOUNDARY PIXELS 

• GROUND TRUTH 

• COLLECTED IN SMALL AREAS AND MAY NOT REPRESENT ENTIRE FIELD 

• TABLE 2 

• FIELDS SHOULD BE DISPERSED TO IMPROVE ACQUISITION PERCENTAGE 

• TABLE 3 

• DIFFERENT PROCEDURES FOR DETERMINING YIELD VARY CONSIDERABLY 


• TABLE 4 
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• VARIATIONS AMONG ACQUISITIONS 

• HAZE 

« SUN ANGLE 
t SOIL AFFECTS 

• GROUND TRUTH SAMPLING 

• LANDSAT IS CORRELATED TO PLANT PROPERTIES THAT INFLUENCE YIELD (TABLE 5) 

• GROUND COVER 

• PLANT HEIGHT 

• STAND QUALITY 

• LANDSAT IS MORE CORRELATED TO YIELD AT SPECIFIC GROWTH STAGES (TABLE 4) 

• JOINTING-HEADING (DAY 107-127 IN EXAMPLE) 



TABLE 1 

SUBSET OF ARS GROUND TRUTH D E P I C T I N G_ C H A NG E S IN YIELD. COMPONENTS 
AFTER OPTIMUM PERIOD FOR SPECTRAL SIGNATURE CORRELATIONS TO YIELD 


Potential Head S1tes/M^ No, Lcaves/M^ 


LOCATION 

ruicnv i o I 

5/17 

ncou 

6/.1 

J 1 LC3/H 

7/23 

I1U« 

5/17 

6/1 

/n 

6/14 

Bu/A 

Seed/Head 

Seed Wt/ 
1000 

10 

157 

278 

264 

583 

948 

809 

20.3 

21.1 

2.43 

1 A 

238 

229 

176 

757 

862 

405 

14.8 

.17.5 

3.01 

6 

378 

370 

346 

1231 

1814 

918 

28.5 

23.0 

2.44 

2 

540 

441 

340 

1836 

1674 

1071 

30.4 

18.0 

3.23 

IB 

819 

730 

408 

2688 

2563 

1958 

3 4. *9 

21.2 

2.66 
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TABLE 2 ! ! 

ARS GROUND TRUTH SHOWING YIELD VARIATION WITHIN FIELDS ! 



YIELD 

(Bu/A) 

LOCATION 

SAMPLE A 

SAMPLE B 

la 

14.8 

27.6 i 

1 

lb 

34.9 

30.9 ‘ ! 

2 

30.4 

21.5 wild oats 
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( • 

24.4 

.25.0 j 

4 

37.9 

42.3 ! 

5 

24.8 

28.6 ' j 

6 

28.5 
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1 7 

32.2 

25.7 
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24.0 . 

29.0 j. ! 

•' ’ 9 
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38.2 

10 

20.3 
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TABLE 3 
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TABLE 

4 - CORRELATION 

(R) OF GREEN 

NUMBER TO: 






YIELO 


NO. 

FIELDS 

ITS 

DATE 

ASCS 

FARMER 

FCIC 

1962 

7531 1 

. 3639 

-.19 


5, 

21 

1963 

7531 1 


. 7934 


8 


1961 

76002 

-.1961 

.5812 

-.1961 

8, 

5, 8 

1962 

76053 

. 2528 

-.2246 


5. 

19 

1962 

7610*7 

-.0159 

.6825 


5. 

21 

1963 

. 76108 


.9303 


8 


1988 

76109 

.7516 

.6169 

. 3993 

23 


1962 

76125 

-.2408 

.4908 


5. 

21 

1964 

76127 

.4473 

.0715 

. 5767 

30 


1988 

76127 

.7256 

. 67C4 

.2170 

23 


196? 

76162 

-.4855 

.0796 


5. 

21 

1964 

76162 

-.0128 

.042 

.2681 



1962 

76197 

-.2984 

.0650 


r>. 

21 

19(3 

76198 


-.4538 


8 


1961 

76200 

. 3729 

-.0508 

.3729 


5. P 

1961 

76234 

- . 5077 \ 

-.5388 

-. 5077 

e. 

5; 8 
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tABLE 


SEGMENT/ 

DATE 

PLANT 

HEIGHT 

GROUND 

COVER 

3 Acqs. 
1964- 76'. 27 
1988-76109, 
16127 

.5416 

.676 

1964 

76127 

.2848 

.2554 

1988 Finney 
76109 

.4502 

• 

.7696 

1988 

76127 

.2548 

.7011 

1988 

76109-76127 

.5367 

.7688 


5 - CORRELATION (R) OF GREEN NUMBER TO : 


GROWTH/ 


GROWTH 

STAGE 

SURFACE 

MOIS- 

TURE 

WEEDS 

FIELD 

OPNS 

YIELD 

DETRAC- 

TANTS 

STANO 

QUALITY 

.4804 

-.1984 

.2131 

-.0775 

.0163 

.575 

.3301 

.3253 


0 

-.1286 

.3526 

0 

.4534 

.0002 

-.2671 

.0177 

.7970 

.0044 

.2632 

.2632 

.0135 

.1617 

.5704 

.464 

-.2057 

.2836 

-.0775 

.2993 

.6302 


(YIELD; 


ASCS 

FARMER 

FCIC 

NO. 

FIELDS 

6194 



76 

4473 

.0715 

.5767 

30 

7516 

.6169 

.3993 

23 

7256 

.6764 

.2170 

23 

6441 

.5621 

.2566 

46 
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SUMMARY 


1 • LANDSAT CAN PROVIDE INFORMATION ON CROP CONDITION 

J 

| • AREAL EXTENT OF MOISTURE STRESS FROM FULL FRAME 

i 

| • SUBJECTIVE RATINGS OF STRESS FROM FULL FRAME CAN BE MADE 

■j . • LANDSAT DIGITAL DATA CAN BE USED TO INDICATE WHEN AGRICULTURAL VEGETATION IS 

}■ UNDERGOING MOISTURE STRESS 

i ' 

I 

| • LANDSAT MAY BE A TOOL TO HELP EXTRAPOLATE PRECIPITATION BETWEEN METEOROLOGICAL 
; STATIONS 

i 

• LANDSAT DATA MAY BE USEFUL IN ESTIMATING SOIL WATER HOLDING CAPACITY 

/ 

• LANDSAT IS SOMEWHAT CORRELATED TO PLANT PROPERTIES THAT INFLUENCE YIELD 

• GROUND COVER 

• PLANT HEIGHT 

• STAND QUALITY 

• LANDSAT APPEARS TO BE CORRELATED TO YIELD AT SPECIFIC GROWTH STAGES ■ 

i 

• ASSESSING YIELD FROM LANOSAT APPEARS FEASIBLE: HOWEVER, MORE RESEARCH IS NEEOEO 



SECTION 6 

ESTIMATING WINTER WHEAT YIELD FROM. 
CROP GROWTH PREDICTED BY LANDSAT 


E. T. Kanemasu 
Kansas State University 


September 27, 1977 
NASA Headquarters 



PREFACE 


The objective of this study is to (1) develop an evapotranspiration 
(ET) model for winter wheat; (2) develop a relationship between Landsat 
data and leaf area index; (3) develop a growth model for winter wheat; 
and (4) develop a yield model using ET and growth models. 

Field data were gathered from commercial fields and plots in 
Riley, Ellsworth, Finney and Thomas counties in Kansas. Data included 
leaf area index, soil moisture, growth stage, and yield. 

Evapotranspiration and growth models required inputs of solar 
radiation, maximum temperature, minimum tempera cure, precipitation, 
and leaf area index. Meteorological data were obtained from National 
Weather Service. Leaf area indices were obtained from Landsat computer 
compatible tapes. Yields were estimated from the ET model; however, 
further testing and evaluation of the yield model are required. 
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1.0 Introduction 


This report summarizes the work completed under NASA Contract 
NAS9-14899. 

Yields are, to a large part, dependent upon solar radiation, temperature, 
and soil moisture. Evapotranspiration and precipitation play = ■ important 
roles la soil moisture. In order to estimate evapotranspiration one requires ‘ 
information as to the vegetative cover. Landsat offers a method of assess- 
ing vegetative cover on repetative basis. Therefore, relatively simple 
weather data supplemented with Landsat estimates of ground cover offer one 
approach to large area yield forecasting. 

2.0 Evapotranspiration (ET) Model 

2 . 1 Model Development 

The daily inputs into the model are solar radiation, maximum-minimum 
temperature, precipitation and leaf area index (LAI). Fig. 1 schematically 
shows the inputs. Potentially, meteorological satellites may be used to 
estimate solar radiation, temperature, and precipitation in areas where 
weather data are not available. Landsat data can be used to estimate LAI. 

The evapotranspiration model described by Kanemasu et al. (1976) 

requires both soil and crop factors to estimate maximum evapotranspiration 

(ET ) and transpiration. ET — the energy- limited ET occurring from a 

max max 

well-watered surface under nonadvective conditions — is given by Priestley 
and Taylor (1972) as 

ET - a[s/(s + y)R 12.1] 

max n 

where a is a constant for a particular crop and climatic situation; y is 
the psychrometer constant (mb/°K) at mean temperature; and is the 24-hr 
net radiation (mm/ day). We evaluated a from lysimetric observations during 
periods of full canopy cover and wet soil surface (a * 1.35). When 
was not measured, ve estimated it from solar radiation, R $ (mm day ^) , using 
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Flov diagram of evapo transpiration (ET) and growth 
models. Potential use of meteorological satellites 
are shown. Winter wheat yields are predicted from 
ET and dry matter production estimates. 










the regression equations: 


R * ,959 R - 3.61 [2.2) 

ns 

and 

R « .926 R - 2.70 [2.3] 

n s 

where [2.2] was developed for growth stages up to jointing and for the 
remainder of the season [2.3). 

Evaporation from the soil surface is limited by energy supplied 

during the constant rate stage; therefore, an energy transmittance term 

t 

(t), based on leaf area index, is required. The daily evaporation rate 
during the constant rate stage can be estimated by 


E - (t/<*)ET 
o max 


[2.4] 


where t * exp(-.737 LAI). Equation [2.4]vas used until ZE * D. Then the 

o 

evaporation was calculated according to the falling rate phase equation 

E f « ct 1/2 - c(t-l) 1/2 [2.5] 

- 1/2 

where c(na day ) depends upon the hydraulic properties of the soil and 
t is days after stage 1 evaporation. The soil factors D and c were 
obtained from lysinetric observations on bare soil or from weight changes 
from large soil-filled containers. 

Transpiration was estimated by equations of the form given by Tanner 
and Jury (1976) and Kanemasu et al. (1976). When the available moisture 
content in the root zone was greater than 352 of field capacity, ve used 


T - o v (1-t)[s/(s + Y>]* n 

crop cover < 502 [2-6] 


and 


T * (a-x)[s/(s +nr))R n 

crop cover > 502 [2.7] 
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where *" 1.56 . 

When the available soil moisture was less than 35Z of the maximum 

available aoisture (Q^). equations [2.6] and [2.7] were multiplied by 
K , given by 

S i 


e a /.35(0 


[ 2 . 8 ] 


Therefore, at 9 less than .35 9 transpiration was linearly reduced as 
a max 

the available water decreased (Fig. 2). The maximum available water content 
of a soil should be determined in the field. 

Soil moisture in the root zone (0-150 cm) was estimated from a water 
balance of evapo transpiration, precipitation, runoff, and drainage. Runoff 
was estimated according to the amount of rainfall (R) and moisture content 
in the surface 30 cm: 


Runoff - 0 R < 2.5 cm [2.9a] 

Runoff - R’ 75 R > 2.5 ca [2.9b] 

where R is the rainfall in inches. The surface 30 cm was allowed to hold 
15 cm of water. Therefore, the rainfall could fill the 30 cm layer to 50Z 
by volume, then the remaining rain must be runoff. The soil profile was 
divided into 5 layers (5, 25, 30, 30, and 60 ca) and each layer was 
allowed to hold 50Z water for two days before draining to field capacity 
(obtained from field measurements). The amount of water drained from the 
5th layer (below 150 cm) was identified as drainage. 

2.2 Procedure 

The evapotranspiration (ET) model was tested on several fields over 
a two year period at Manhattan, Kansas. Daily estimates by the model were 
compared with lysimetric observations. Leaf area index (LAI) was measured 
by optical planioeter and/or leaf length and width calculations. Soil 
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Fig. 2. Water stress factor (Ks) as a function of available 
water in the root zone. Ks linearly declines at 
3SZ available water. 
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moisture estimates by the model compared favorably with neutron attenuation 
and gravimetric estimates* 

LAI obtained from ground measurements are extremely tedious. Landsat 

i 

data were used in the ET model by estimating LAI* Multiple regression 
equation was developed from Landsat coverage of Kansas sites (Colby, 

Ellsworth and Manhattan, Table I)* Shown in Fig. 3 is the comparison of 
Landsat-predicted LAI with observed LAI. Figs. 4 and 5 show the season 
trends in observed and Landsat-predicted LAI. When Landsat predicted LAI 
curves were used in the ET model instead of observed LAI, seasonal ET 
estimates by Landsat were usually within 3.0 cm of the ET estimates from 
observed LAI measurements. 

3.0 Soil Moisture Estimates from ET Model 
For the 1975-76 winter wheat growing season, we obtained sample 
statistics for 22 sample segments in five Great Plains states (Kansas, Texas, 
Oklahoma, Nebraska, and Colorado). Analyst interpreters selected several 
wheat fields in each segment (4 to 20 fields). Landsat data were analyzed 
for each useable overpass date on all fields. For each date, leaf area 
index was estimated for each field and then averaged to obtain an average LAI 
for the segment (Figs. 6 and 7). The ET model was run on each segment and 
estimates of soil water depletion (higher percent depletions are drier) 
throughout the growing season are predicted (Figs. 8, 9, 10). 


4.0 Yield Estimates from ET Model 

A yield model was developed from small plot yields and the output from 
the ET model. 


Yield (metric tons/ha) * 


(Z(T/ET max )] 


0.104 

2 


0.192[£(T/ET max )]J' 172 * 

[E(T/ET ))°" 646 
1 max 3 


[4.1] 
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Table 1. Computer compatible tapes 
used la data analysis. 

TAPES USED IN 

AREA DATE TAPE I.D. # 

Colby 8/20/75 5123-16310 
8/29/75 2219-16442 

9/07/75 5141-16300 

9/25/75 5159-16285 

10/22/75 2273-16440 

1/11/76 5267-16221 

2/25/76 2399-16421 

4/01/ 75 2435-16410 

4/10/76 5357-16161 

6/02/76 5410-16065 
6/03/76 5411-16123 

6/12/76 2507-16391 

6/20/76 5428-16053 

6/30/76 2525-16384 

7/09/76 5447-16095 

9/10/76 2597-16364 

10/16/76 2633-16353 

11/21/76 2669-16341 

AREA DATE TAPE I.D. # 

Manhattan 10/20/73 1454-16374 

3/31/74 1616-16344 

4/18/74 1634-16341 

5/24/74 1670-16331 

6/29/74 1706-16320 

7/17/74 1724-16313 

8/04/74 1742-16305 

9/09/74 1778-16293 

10/15/74 1814-16283 

11/20/74 1850-16272 

12/07/74 1867-16205 

3/25/75 1975-16161 

4/12/75 1993-16152 

4/30/75 5011-16142 

5/18/75 5029-16133 

6/06/75 5048-16181 

6/24/75 5066-16171 

8/16/75 5119-16082 

11/15/75 5210-16083 

12/03/75 5228-16073 

4/16/76 2540-16232 

5/04/7-: 2468-16225 

6/09/76 2504-16220 

6/17/76 5425-15483 

7/06/76 5444-15525 

9/06/76 2593-16135 

9/24/76 2611-16131 

10/12/76 2629-16124 

10/13/76 2630-16182 


from Landsat multispectral scanner 


DATA ANALYSIS 

, AREA DATE TAPE I.D. # 

Ellsworth 9/23/75 5157-16173 

10/02/75 2253-16324 

10/11/75 5175-16163 

10/20/75 2271-16323 

10/29/75 5193-16152 

11/07/75 2289-16322 

11/16/75 5211-16141 

1/18/76 2361-16313 

3/12/76 2415-16301 

3/21/76 5337-16061 

3/30/76 2433-16294 

6/01/76 5409-16011 

6/10/76 2505-16274 

7/07/76 5445-15583 

10/14/76 2631-16240 

11/01/76 2649-16233 

11/19/76 2667-16224 

12/25/76 2703-16211 


AREA DATE TAPE I.D. # 

Manhattan 10/31/76 2648-16174 

11/17/76 2665-16112 

11/18/76 2666-16170 

12/24/76 2702-16153 
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OBSERVED LAI 








SOU. WATER DEPLETION (X) 


1975-76 

WINTER WHEAT 
GA*0€\‘ CO, NEBRASKA 



DATE 


7. Seasonal trends in Landsat-predicted leaf 'area 
index (LAI) for sample segment in Garden County 
Nebraska , 1975-1976 . 



DATE 


. 8. Seasonal trends in soil water depletion in Grand 
County, Kansas, 1975-1976. 





Fig. 9. Seasonal trends in soil water depletion in Kearney 
County, Kansas, 1975-1976. 



Fig. 10. Seasonal trends in soil water depletion in Stevens 
County, Kansas, 1975-1976. 
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where the subscripts l f 2, and 3 are the respective growth stage intervals: 
emergence to joxnting, jointing to heading, and heading to soft dough; T 
is the daily transpiration rate; ET m is the energy-limiting evapo- 
transpiration rate. Therefore, the yield model can be used on any field 
where the ET model can be applied. 

Eleven wheat fields at Bushland, Texas presented an independent data 
set. Landsat and yield data were available (personnal communication with 
Dr. Clif Harlan, Texas A & M) . The ET model was run using meteorological 
data and Landsat-predicted LAI. Yields were predic ed from [4-1] and 
compared with observed yields (Fig. 11). 

The soil moisture study over the 5 Great Plains states offered 

. 

another data set; however fields for individual fields were not measured. 
County yields were available from the Statistical Reporting Service (SRS) . 
in addition, Feyerherm's KSU winter wheat model was run on the same data 
assuming a management and productivity (MAP) factor of 1 and summer fallow 
conditions. The root mean square error (RMSE) between the county yield 
and the ET yield model (eq. [A.l]) was 2.0 bu/acre while the RMSE between 
Feyerherm's yield model and the ET yield model was 1.5 bu/acre. 


5.0 Growth Model 

As shown in Fig. 1, the growth model uses the identical inputs as 
the ET model — solar radiation, max-min temperature, precipitation, 

LAI. The major assumption in the growth model is that light and P~il 
moisture are the primary limiting factors in plant growth. Other factors 
such as fertility, Fest and disease influence growth and are reflected In 
the LAI term. 

Photosynthesis is estimated from the amount of light that the canopy 
intercepts which is dependent upon the solar radiation and LAI. Soil 
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moisture decreases photosynthesis during high water depletion periods. 
Respiration is dependent upon 1AI and temperature. The difference between 
photosynthesis and respiration is net photosynthesis which is the rate of 
dry matter production . The growth model simulated dry matter 

production on commercial fields in western, central and eastern Kansas 
using measured LAI. Fig. 15 shows the agreement in dry matter production 
estimated by the growth model using Landsat-predicted LAI and observed 
LAI. _ - - 
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DRY MATTER (MEASURED LAI) 


. Comparison of dry matter estimated by the growth 
model using measured leaf area index (LAI) and 
Landsat-predicted LAI. 
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SECTION 7 


A STUDY OF LASDSAT DATA 
AND 

EARTHSAT SPRING WHEAT DATA 
FOR 

YIELD DETERMINATION 


Executive Summary 
Adapted froa Contract KAS5-22950 
July 1976 Goddard Space Flight Center 


Prepared For 
THE 

LANDSAT CROP CONDITION 
AND 

YIELD BRIEFING 
September 27, 1977 


By 

Earl S. Merritt 
Vice President 
Earth Satellite Corporation 
7222 47th Street 
Washington, D. C. 20015 



EXECUTIVE SUMMARY - STATEMENT OF EARTHSAT INTERACTIVE YIELD ESTIMATE 

CONCEPT 

The EarthSat Yield System has been developed as a modem alternative 
to the traditional weather regression approaches to crop yield estimates. 
The "System" has, furthermore, been designed from inception to permit 
the interactive use of yield-related information derived from remote 
sensor systems, either aircraft or satellite. 

The EarthSat System is largely computerized. It operates on a 
global ly-appliceble two-level (25n.ra. and 12.5n.n.) geobased gnd-cell 
structure. The "System" processes meteorological data from first order 
ground meteorological observation stations and from meteorological 
satellites in order to define a dense network of real and synthesized 
plant weather information. In the 1975 upper Great Plain tests, weather 
station data and meteorological satellite data were entered into the 
"System" at six hourly intervals. 

The objective of the basic diagnostic activities in the "SystemT is 
to define the weather influencing plant growth with sufficient detail 
that simulation models which describe plant growth, and define soil 
moisture profiles can be accurately operated. The goal of all system 
diagnostic activities is to define the spatial variations in plant yield 
clearly enough that such descriptions can be locally verified with 
either ground-based observer transects or by remote sensing techniques. 

The "System" differs from traditional approaches in that the re- 
sulting synthesized and real weather diagnostic crid allows application 
of quasi-physiologically and fully physiolcgically-based plant yield 
models. These models either describe or infer plant processes, i.e. , 
photosynthesis, gas exchange, dry matter accumulation and translocation. 
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water stress, etc., accurately enough to permit a very accurate and 
highly plant descriptive diagnosis. 

The plant process descriptions over a 12.Sn.in. geobased cell 
structure for the 1975 spring wheat crop season has been utilized to 

9 

develop a functional relationship between LANOSAT observables and the 
stress factors described by the "System." This functional relationship, 
which includes a component of the short term plant stress as well as the 
long term stress history, has been utilized to enhance the spring wheat 
yield forecasts produced by the simulation model. These enhanced yield 
estimates were prepared after complete analysis of all LANOSAT frames 
taken between 15 Xay and 15 September over the four state upper Great 
Plains region. 

The LA.NDSAT analyses were undertaken using a previously defined 
interpretive key which permitted a description of low, moderate or high 
stressed areas with an approximate 65 percent accuracy and lew and high 
stress area with an accuracy of approximately 90 percent. The LANOSAT 
analyses were then coded for entry into the computerized geobase at a 
resolution of approximately 12.5n.m. Once entered into the data base 
they were readily available for interactive uses with the existing 
EarthSat “System'' 1 simulated data base. 

The results achieved by the LANOSAT Interactive EarthSat System 

show definite promise. For example, at the four state aggregate level 

the error of yield estimate was reduced from an already reasonable 2.35 

to 0.79%. At the state level the average error of approximately 55 
* 

was lowered to an average error of approximately 35. Similar improve- 
ments were generally noted at the crop reporting district (CRO) level. 
The region-wide errors produced by NOAA's traditional regression models 
for the same area and time were 6.35. Table £-1 presents these com- 
parisons. 
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ORIGINAL PAT IK IS 


EARTHSA1 YIELD ESTATES 



BASED 0.1 DATA THRU AUGUST >3 1975 

m xs:m 

mt Mk -aifsirH 

WAA 

CCLA 

8/25/75 

ESTi^Tc 

m 

3/10/76 

M* 

HSITASt 

23.7 

23.01. 

21.39 _ 

_ -2Z.6 _ 

25.6 

a.c. 

2S.0 

29.31 

21.96 

29.S 

27.2 

a.E. 

23.3 

22.63 

21.16 

26.7 

2S.0 

c. 

— 

21.65 

21.99 

— 

25.6 

s.c. 

— 

— 

— 

— 

20.1 

S.E. 

— 

21.37 

19.59 

— 

23.7 

SSHR RAaDTA 

17.0 

33.31 

_ 17.35 

70.7 

13.0 

ii. y. 

13.7 

19.31 

13.33 

20.3 

17.9 

H.C. 

15.9 

17.32 

15.23 

19.5 

19.0 

a.E. 

17.3 

21.19 

17.56 

13.1 

19.9 

w.c. 

— 

. 


— 

13.7 

c. 

— 

19.C3 

16.11 

— 

15.9 

E.C. 

— 

19.61 

13.72 

. 

17.9 

s.y. 

— 



_ . _ 

— 

15.0 

s.c. 

— 


^ , 

— 

15.3 

S.E. 

— 

18.09 

16.99 

— 

17.9 

RORTH PASTA 

27. X 

2S.92 

2G.95 

73 A 

25.9 

H.V. 

25.6 

29.39 

25.9S 

21.1 

29.9 

H.C. 

26.2 

29.79 

29.99 

20.2 

29.9 

II. E. 

29.7 

29.53 

29.03 

29.6 

31.0 

H.C. 

— 

29.99 

26.36 

— 

29.6 

C. 

26.7 

26.22 

25.82 

22.0 

25.5 

E.C. 

33.7 

29.91 

29.33 

30.1 

23.3 

s.y. 

— 

23.13 

22.21 

— 

22.9 

S.C. 

— 

23.26 

23.92 

— 

21.2 

S.E. 

23.3 

27.33 

26.83 

16.3 

22.8 

HIIHIESQIA. 

33.7 

33.93 

_33J1L 

30.2 

.31 JL 

a.y. 

36.2 

39.35 

39.50 

33.5 

33.7 

H.C. 

— 



— 

26.5 

H.C. 

29.5 

32.99 

30.65 

ZlA 

27.0 

c. 


— 

m 

— 

31.0 

s.w. 

— 

29.32 

27.71 

— 

33.0 

s.c. 

■ 

— 



— 

31.0 

Foua State 
ESTIMATE 

26.0 

25.6 

25.2 • 

23.S 

25.9 


• LARDSAT CORRECTED STRESS HISTORY iiODEL 

•' IilTERACT lVS SOIL fwlSTURC AiJD UWOSAT CORRECTED STRE SS H ISTORY liODEL 

TABLE t-TH 
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The LAfiOSAT analytical technique has been applied to winter wheat 
areas of Kansas and surrounding states in 1975 and 1976. In this period 
concern over a new "dust bowl" in northwest Oklahoma, southwest Kansas, 
northwest Texas and eastern Colorado was high in the late winter and 
early spring of 1976 since poor germination had been observed over much 
of the area. The resultant LANOSAT analysis accomplished in October 
1975 and March through )May 1976 indicated that, for the state of Kansas 
the fears of a “dust bowl" were only justified over southwest CRD of 
Kansas where extensive abandonment of dry land winter wheat fields 
occurred. All other areas of Kansas were reasonably good but they were 
be' 1 their record 1975 yields. Total production was down nearly 71.5 
million bushels over 1975. 

The full EarthSat Interactive "System" was not operated over the 
winter wheat region. However, the system models were operated from 
planting to 1 April at selected ground observation points. These sample 
runs appear to confirm the applicability of the "System" diagnostic and 
predictive element in the winter wheat areas. Selected point average yield 
estimates are Dodge City 23bu/ A, Topeka 24bu/A, Amarillo 14bu/ A. These 
yield estimates are based on the use of a Technology Acceptance maximum 
yield value derived from the past .4 years of Kansas yield history and 
plant stress coefficients developed in the spring wheat region states. 

The EarthSat Yield System concept has shown considerable promise in 
the spring wheat test in 1975. The use of LANDSAT interpretation gen- 
erally appears to improve the "System" yield estimate. The application 
# 

of all types of data in a common coordinate system is a very powerful 
concept. The combination of this concept with a highly disaggregated 
plant environment diagnostic and plant yield simulation (process) models 

1/ Includes both dryland and irrigated area yields. 
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offer additional improvements in the future. It is anticipated that the 
greatest benefits from the EarthSat System will accrue to yield estimates 
made in anomalous years and in regions where the meteorological observing 
network is less dense than in the United States. 

EarthSat CROPCAST™ System, a commercial crop forecasting venture, employs 
some aspects of the System studied in 1975 and 1976. CR0PCAS7 is now in opera- 
tion over Canada and the United States for corn, soybeans, wheat and cotton. 
Results to date are encouraging, e.g., comparisons of CROPCAST's forecasts 
of the USDA monthly (SRS) Crop Production Reports, issued approximately 4 weeks 
and two weeks prior to the USDA report, show the following accuracies: 


All Crops 

97% 

Corn 

98% 

Soybeans 

97% 

Winter Wheat 

99% 

Spring Wheat 

93% 


The end-of-year comparisons are a few months away, but similar accuracies are 
expected. 

CROPCAST is now available over the South American soybean and wheat areas 
in Brazil, and Argentina. Monitoring of Winter conditions is underway over 
several wheat growing regions. 

CROPCAST has been designed to use Landsat when it is available in a timely 
manner. The future plans for 48 to 95 hour turn-arounds are very exciting. 
CROPCAST will continue to use Landsat in a confirmatory and interactive manner, 
rather than as a primary data source. 
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SECTION JL- 

SPECTRAL INDICATORS OF CROP DEVELOPMENT AND LEAF AREA INDEX 

FROM LANDSAT DATA 

C. L. WIEGAND, H. W. GAUSMAN, A. J. RICHARDSON, 

A. H. GERBERMANN, J. H. EVERETT, AND R. W. LEAMER 

Abstract 

Spectral indices such as the transformed vegetation index (TVI), 
the green number (GVI), and the perpendicular vegetation index (PVI) 
are significantly correlated with leaf area index (LAI), and green 
biomass (BIOM) during the crop development and grain filling stages. 

They also respond to growing conditions as LAI and BIOM do. Two of 
them take soil background into account, hence also help remove its vari- 
ations in MSS data. By so doing, they offer the possibility of calibrat 
ing crops spectrally across years, thereby minimizing ground truth 
requirements and increasing the value of the indices where ground truth 
is unavailable. In addition, they and their temporal trajectories may 
be helpful in improving training sample selection, signature extension, 
and in classification procedures. 

The evidence indicates that the vegetation indices can be used to 
estimate LAI needed for the evapotranspiration and photosynthesis 
subroutines in crop productivity models. Thus they can be used to help 
implement the models over large areas by either (a) providing input data 
for the models, or (b) feedback da^a to check on, and retrack the models 
if necessary. 
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LAND SAT FOLLOW-ON FINAL REPORT 
CONCLUSIONS 

1. POP, PC, MID LAI ARE THE PLANT PARAMETERS HOST CONSISTENTLY 
RELATED TO LAND SAT MSS DIGITAL COUNTS (DC), 

• LAI CAN BE ESTIMATED SPECTRALLY. 

• LINEAR COMBINATIONS OF THE OTHER PLANT PARAMETERS (POP, PC, 

PH) ACCOUNT FOR 67 TO 90% (R J X100) OF THE VARIATION IN LAI 
AND FROM 69 TO 89% OF THE VARIATION IN GRAIN YIELD. 

2. LANDSAT SPECTRAL INDICATORS, SUCH AS FVI, RELATE TO GRAIN YIELDS 
OF SORGHUM FOR ABOUT A 60-DAY PERIOD — FROM GROWING POINT DIF- 
DERENTIATION (GPD) TO HALFWAY BETWEEN 1/2 BLOOM (HB) AND PHY- 
SIOLOGICAL MATURITY (PH) OF THE GRAIN. 

3. OPTIMAL WAVELENGTHS FOR DETECTING CERTAIN STRESSES HAVE BEEN 
DETERMINED. 

*». FORAGE PRODUCTION DIFFERENCES OF GRASSY RANGELANDS CAN BE HAPPED. 
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"SPECTRAL INDICATORS OF SORGHUM DEVELOPMENT AND THEIR 
IMPLICATIONS FOR GROWTH MODELING" 

CONCLUSIONS 

(f TEMPLE, TX 1976 SORGHUM DATA) 

VEGETATION INDICES DERIVED FROM LAND SAT DATA ARE RESPONSIVE TO 
GROWING CONDITIONS THAT AFFECT LAI AND BIOMASS. 

TV I, PVI, and GVI are about equally useful for monitoring 
seasonal crop development and vegetation density. 

THE HIGH CORRELATIONS OBTAINED BETWEEN LANDSAT VEGETATION INDICES 
AND PLANT GROWTH MEASUREMENTS INDICATE THEY CAN BE USED OVER LARGE 
AREAS, EITHER AS 

a) INPUT DATA FOR PLANT GROWTH SIMULATION, OR 

b) FEEDBACK DATA TO CHECK ON, AND RETRACK GROWTH SIMULATION 



3. SPECTRAL VEGETATION INDICES CAN BE CALCULATED FOR AS MANY 
PIXELS, OR FIELDS, AS ARE OF INTEREST IN A GEOGRAPHICAL 
AREA. THUS, PLANT GROWTH MODELS CAN BE EXTENDED TO LARGE 

• AREAS YET BE AIDED BY SPECIFIC FEEDBACK ON ACTUAL GROWING 
CONDITIONS IN INDIVIDUAL FIELDS. 

4. THE IMPROVED ESTIMATES OF LEAF AREA INDEX AND BIOMASS THAT 
RESULTED FROM INCLUSION OF WEATHER DATA IN COMBINATION WITH 
VEGETATION INDICES IN ESTIMATING EQUATIONS INDICATE THAT GROWTH 
SIMULATION MODELS THAT MIMIC PLANT RESPONSE 70 SOIL AND AERIAL 
ENVIRONMENTS WILL IMPROVE YIELD ESTIMATES* OVER THOSE ARRIVED 
AT FROM SPECTRAL DATA ALONE. 


* FARMER-REPORTED YIELDS ARE SUSPECT! (DISAGREE WITH BOTH GROUND 
SAMPLE DATA AND SPECTRAL INDICATORS.) 
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Sorghum Plant Growth Measurements 
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** Significant at the 0.01 probability level. 
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"LEAF AREA INDEX ESTIMATES FOR WHEAT FROM 


LAND SAT SPECTRAL DATA" 

(Wiegand , Richardson, Kanenasu ) 

CONCLUSIONS 

1. FVI = *VCRgg5 - RpS) 2 + (Rgg7 - Rp7) 2 

YIELDED EQUAL OR BETTER CORRELATION WITH GROUND -MEASURED 
LAI THAN DID 

LAI = a Q - a^(HSS 4/5) - a 2 (MSS 4/6) + a 3 (MSS 4/2x7) + a 4 (MSS 5/6) 
- a 5 (MSS 5/ (2x7)) + a 6 C(MSS 4/5) - (HSS 4/(2x7)] MSS(4/5) 

2. SPECTRAL VEGETATION INDICES SUCH AS PVI ARE APPLICABLE TO 

3. APPEARS POSSIBLE TO CALIBRATE WHEAT LAI IN TEi A S OF PVI AND REDUCE 
GROUND TRUTKING TO SPOT CHECKS. 
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